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Abstract- In a regulated environment, unit commitment refers to 
optimizing generation resources to satisfy load demand at least 
cost. But in a deregulated power market, generation companies 
have no obligation to satisfy customers' demands. Instead, they 
seek their maximum profits. This paper proposed a price-based 
unit commitment algorithm considering uncertainties for 
generation companies to maximize their profits. Due to 
uncertainties existing in energy prices, a fuzzy regression model 
is proposed to handle these uncertainties. Based on the fuzzy 
energy price model, an algorithm considering uncertainties of 
energy prices solved by cooperation of the greedy algorithm and 
quadratic program is presented in this paper. The New York 
electrical market data are used to verify this algorithm, and test 
results of the 5-unit system and 36-unit system are presented in 
this paper. 

Keywords- Deregulated Power Market; Price-based Unit 
Commitment; Fuzzy Regression Model; Greedy Algorithm; 
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I. INTRODUCTION 

For the vertically regulated power industry, unit 
commitment (UC) refers to optimizing generation resources to 
satisfy load demand at least cost. But in the deregulated power 
market, generation companies (GENCOs) run UC not for 
minimizing operation costs but for maximizing profits for their 
own. Satisfying load is no longer an obligation for the 
GENCOs, and security would be unbundled for energy and 
worked as an ancillary service. The GENCO can only consider 
a schedule that obtains the maximum profits neglecting the 
load demand. 

Many methods [1-9] have been proposed for the traditional 
cost-minimization unit commitment. In the competitive 
environment, UC plays a new role to maximize GENCOs' 
profits. At the stance of GENCOs, new formulations have been 
proposed to solve this problem. In [10], the authors formulate a 
profited-based unit commitment considering power and reserve 
generating. A hybrid method between Lagrange Relaxation and 
Evolutionary Programming is applied to solve this new UC 
problem. In [11], the authors describe a genetic algorithm 
solution to the new profit-based unit commitment problem. 

In [12], the authors formulate the price-based unit 
commitment (PBUC) problem based on the mixed integer 
programming (MIP). The PBUC solution by utilizing MIP is 
compared with that of Lagrangian relaxation method. Test 
results show the efficiency and advantages of the MIP methods 
for solving PBUC. In [13], a security-constrained price-based 
unit commitment (SPUC) algorithm underlining the 
coordination technique which is implemented in two stages is 
proposed. This paper describes a coordination process between 
GENCOs and the ISO for congestion management and 
reducing the risk of failure to supply loads. At first, GENCOs 
apply priced-based unit commitment (without transmission 
security constraints), schedule their generating units and submit 



their energy and capacity bids to the ISO for maximizing their 
revenue. The ISO obtains transmission information as well 
from TRANSCOs via OASIS. Using the two sets of generation 
and transmission data, the ISO executes inter-zonal and intra- 
zonal congestion management and contingency analysis to 
minimize line flow violations and the risk supplying loads. In 
[14], a methodology based on dynamic programming and 
enumeration has been proposed to solve the price based unit 
commitment problem. The technique has further been extended 
to optimize the problem for a group of generators. 

Because energy prices are unknown before bidding, 
uncertainties of energy prices are important factors to build the 
new UC problem. Therefore, this paper presents a price-based 
UC problem considering uncertainties of energy prices at the 
stance of GENCOs. A fuzzy regression model is used to take 
uncertainties exiting in energy prices into account. The New 
York electrical market data are adopted to verify this model. 
This problem is solved by the greedy algorithm and quadratic 
program using the concept of decommitment [9]. For GENCOs, 
the unit may be turned off in some on-line periods only if doing 
so is cost-saving and would not cause infeasibility. Given a 
feasible unit commitment, this algorithm uses the greedy 
algorithm to determine an optimal strategy for decommitting 
overcommitted units. This paper is organized as follows. The 
energy price model is described in Section II. Section III will 
introduce the price-based unit commitment formulation. The 
greedy algorithm and quadratic program used to solve this 
problem are specified in Section IV. Finally, test results of the 
5-unit and 36-unit systems are given in Section V. 

II. THE ENERGY PRICE MODEL 

Due to uncertainties of the energy prices, forecasts of 
energy prices are essential information to the GENCOs. In this 
paper, a fuzzy regression model [15] is used to forecast the 
range of energy prices, and a method to build the energy price 
scenarios of the next day's twenty-four hours is proposed 
according to the forecast ranges of energy prices. 

The most important factor to the energy price is the load 
demand. One of the most popular methods for representing the 
relation between these two sets of data is regression analysis. In 
regression analysis, the relation between parameters x and y is 
represented by the following polynomial equation: 

y = a +a 1 x + a 2 x 2 + + a m x m (1) 

By choosing the order m we can represent nonlinear 
relations. Parameters are determined so that the distance (or 
error) between an observed data point and its corresponding 
point on the polynomial will be minimal. 

Fig. 1 shows the demand-price data of the New York 
electrical market from January 2002 to May 2005. When the 



IJ ES Vol.2 No. 2 2012 PP. 51-58 www.ijesci.org ©World Academic Publishing 

~ 51- 



International lournal of Energy Science 

demand is low, the relation between the price and demand is 
much like linearity. However, as the demand increases, the 
trend pushes up the price. Therefore, a linear regression model 
is applied to low demand data, while higher order curves are 
more suitable in representing high demand data. 

Price (USt) 




5000 10000 15000 20000 25000 30000 35000 
Demand (MWh) 



Fig. 1 The demand-price data of the New York electrical market 

The TSK-fuzzy model was developed by Takagi, Sugeno, 
and Kang [16-17] to represent the nonlinear relation of multiple 
inputs and output data, based on the format of fuzzy reasoning. 
In this paper, we divide the data clusters of demand into two 
overlapping regions-low demand and high demand-such as 
shown by the membership functions of fuzzy sets low and high 
in Fig. 2. 

tot 




Fig. 2 Fuzzy sets to represent low and high demand 

The low and high data clusters are supposed to overlap 
between D L and D H . Linear and quadratic regression curves are 
applied to low and high demand data respectively. Such a fuzzy 
model may be represented by the following rules: 

IF x is low THEN yj =ag' + aj''x n) 
IF x is high THEN y 2 = a{, 2 > + aj 2) x + a ( 2 2) x 2 

where in this case x represents the demand and y-„ i=l, 2 stand 
for electricity prices. 

The composite nonlinear input-output relation is obtained 

by 

y = Eft-yj (3) 
I Hi 

where is the degree of membership of the demand data 
belonging to either the low (i=l) or high (i=2) fuzzy set. 

When x is smaller than D L , linear regression is solely 
applied. When x is greater than D H , only quadratic regression is 
applied. When x falls in between, both equations are employed 
with the continuously varying degree of weight . 

A. Fuzzy Regression Model 
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In the fuzzy regression model, the parameters in (1) are 
replaced with fuzzy numbers as shown in (4) to cover a wide 
range of data. 

Y = A +A 1 x + A 2 x 2 + --- + A m x m (4) 

The parameters Ao, Ai,....are so determined that the 
observed data are encompassed by the fuzzy regression model, 
and the resultant left-hand variable Y is also a fuzzy number, 
which has a region of data covered in a varying degree of 
possibility. Fig. 3 shows the triangular fuzzy set representing 
the fuzzy number Aj with three crisp parameters, namely 

Cj , Cj . 



1 








Fig. 3 Triangular fuzzy set to represent regression parameters 

Here, a; is the most likely value of the regression parameter, 

whereas c^ and cf are possible maximum spread from a; to the 
higher and lower values of the parameter, respectively. We use 

the expression A; = (a^c* ,Cj r ^to represent such a triangular 

fuzzy number hereafter. In the modeling process, the mean 
value a ; of the fuzzy number is most simply determined by 
conventional regression, but, naturally, it can also be obtained 
by other methods such as neural networks. 

To determine the remaining parameters for the fuzzy 
numbers (cf and cj"), a linear programming model is used to 
fit the model to the given data within as small an area as 
possible. The optimization process is formulated as follows: 

Minimize 

i[; ( ;+cfx(t) + - + c;x(t) m +Co+c[x(t) + -+c^x(t) m } 
t=l 

S.T. 

a +a,x(l) + a,x(l) 2 + • ■ ■ + a m x(l) m + cj + cfxfl) + cjx(l) 2 + ■ ■ • + C ;,x(l) m > y(l) 

: " (5) 

a +a 1 x(n) + a 2 x(n) 2 +---+a m x(n) m +cj +c ] f x(n) + cjx(n) 2 +---+cj,x(n) m > y(n) 
a + a 1 x(l) + a 2 x(l) 2 + --- + a m x(l) m ^ 

a +a 1 x(n) + a 2 x(n) 2 + --- + a m x(n) m -cjj -cj"x(n)-c 2 x(n) 2 c m x(n) m < y(n) 

Co.Co.C^C] ,-,c m ,C m >0 

where n is the total number of observed data points. 

In order to narrow the range of the fuzzy regression model, 
some extreme market data that have little chance happening are 
eliminated. The modified price-demand data are shown in Fig. 
4. Applying these data to the above model, parameters for the 
low and high demand are described in Table I and II. 
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Price (USf) 




5000 10000 15000 20000 25000 30000 35000 
Demand (MWh) 



Fig. 4 The modified demand-price data of the New York electrical market 



TABLE I PARAMETERS FOR LOW DEMAND MODEL 





1 st Order 


Constant 


a 


3.27484xl(T 3 


-13.23068 


c + 


3.3579xl0" 3 





c" 


1.7781xl0" 3 






TABLE II PARAMETERS FOR HIGH DEMAND MODEL 





2 nd Order 


1 st Order 


Constant 


a 


1.57089xl0 7 


-1.708595xl0" 3 


-1.70383x10 3 


c + 


7.32589x10 s 








c" 


4.97407x10 s 









Parameters for Table I and II are built by the demand-price 
data from January 2002 to May 2005. Using these parameters 
and the demand data of June 1 to 7, energy prices of June, 1 to 
7 can be estimated as shown in Fig. 5. 



Actual — ■ — Crisp — ^ Up (a = 0) Down (a = 0) 




06/01 06/02 06/03 06/04 06/05 06/06 06/07 
Fig. 5 Energy prices estimated by fuzzy regression model 



Even after reasonable data modification, the resultant range 
of fuzzy model may appear too broad because the optimization 
process introduced above considers the entire possible range of 
the observed data. One idea to narrow down the focus in an 
integral fashion is to apply a - cuts of the fuzzy numbers in (4) 
to modify their shape. Using a single parameter a(0<a<l), 
the base of the triangular fuzzy numbers are modified as the 
equation (6). At a = 0, we have the original fuzzy regression 
model. As the value of a increases toward 1, the range 
encompassed in the fuzzy model narrows, and eventually it 
shrinks to the convention regression model at a = 1 . For 
a > 0, some of the data points outside the range of the fuzzy 
model are neglected, and the number of missed data will 
increase as a approaches to 1 . Therefore, the value of a can be 
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interpreted as indicating the willingness of the decision-maker 
to accept the risk of missing data points to obtain a narrow 
range. 

Aj ={& i ,c\ x (1 - a), cj~ x (1 - ot)y (6) 

B. Energy Price Scenarios 

Due to uncertainties of energy prices, different energy price 
scenarios are considered in the price-based unit commitment 
model. The energy price scenarios are established according to 
the estimated energy price range. They are built by choosing 
a value from 0.5 to 1.0 with the step of 0.1. Except 
a = 1 .0 with one curve, each a value corresponds to two 
curves. Therefore, 1 1 energy price scenarios are created. Each 
energy price scenario is given a probability proportional to its 
a value by the following equation. 

P (7) 

s s 

i=l 

where 

P s is the probability of the energy price scenario s, 
a s is the a value of the energy price scenario s, 
S is the number of energy price scenarios. 

Fig. 6 shows the 11 energy price scenarios of June 1, 2005 
according to the fuzzy regression model. The probability for 
each curve is calculated by equation (7). These energy price 
scenarios will be used by the proposed price-based unit 
commitment model. 

Price(SUS) 





— ♦ P s =0.125 

— ■ — P s =0.1 125 

P s =0.1125 
t P s =0.1 






9 P s =0.0875 

1 P s = 0.0875 

P s = 0.075 

■ P s = 0.075 

P s - 0.062: 
P s =0.0625 




2 4 6 8 10 12 14 16 18 20 22 
Time 



Fig. 6 Possible energy prices scenarios (2005/06/01) 
III. THE PRICE-BASED UNIT COMMITMENT MODEL 

This new price-based UC problem is different from the 
traditional cost-minimization UC. The goal of the price-based 
unit commitment is to maximize the profit (i.e., revenue minus 
cost) of GENCOs subject to all constraints. In the deregulated 
environment, energy prices are unknown when GENCOs 
propose their bids. Therefore, this new algorithm should 
consider uncertainties existing in energy prices and the 
objective function is to maximize profit instead of cost 
minimization. Besides, the new UC problem doesn't consider 
satisfying load demand as its necessary constraint. The 
objective function and constraints of this model are described 
below. The objective function of this model is first depicted as 
follows. 

max £P S *(iZRV(s,i,t)-TC(i,t)) (8) 

s=l i=lt=l 

where 

RV is the revenue of GENCO 
TC is the production cost 
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S is the number of energy price scenarios, and P s is the 
probability for each scenario. 

RV(s, i, t) = PE(s, t) * P(i, t) * I(i, t) (9) 

TC(i,t) = {C i (P(i,t)) + St(i,t)}*I(i,t) (10) 

where 

PE(s,t) is the energy price at hour t for scenario s, 
P(i,t) is energy generation of generator i at time t, 
St(i,t) is the start up cost, 
I(i,t) is the unit status of unit i at time t, 

C ; (P(i, t)) = a(i) + b(i) * P(i, t) + c(i) * P(i, t) 2 
is the production cost of unit i at time t, 

S is the number of energy price scenarios, and P s is the 
probability for each scenario. 

From equation (9), it can be seen that the revenue is from 
the sales of energy. Profit is defined as the revenue minus 
production costs. Constraints for the new UC problem are also 
different from those of the traditional one. These constraints 
are described below. 



(a) System energy limits 

£p(i,t)*I(i,t)<D t , t=l,-T 
i=l 

where 

D t is forecasted demand at hour t. 

(b) Unit generation constraints 



1 grain 



(i)<P(i,t)*I(i,t)<P gmax (i) 



(ID 



(12) 



where 



1 gmin 



(i) is minimum generation limit of generator i, 



Pgmax (i) I s maximum generation limit of generator i. 
(c) Unit minimum on/off durations 

t+T on (i)-l 

I I(i,t)>a it *T on (i) t = l,...,T-T on (i) + l 
t 

T 

£I(i,t)>a it *(T-t + l) t = T-T on (i) + 2,...,T 
t 

t+ T off (i)-l 

X [1 - I(i, t)] > P it iT off (i) t = l,...,T-T off (i) + l 



I[l-I(i,t)]>p it *(T-t + l) t = T-T off (i) + 2,...,T 



(13) 
(14) 
(15) 
(16) 



where 



a it is a binary variable that is equal to 1 if unit i is started 
up at hour t and is otherwise, 

(3 it is a binary variable that is equal to 1 if unit i is shut 
down at hour t and is otherwise, 

T is the planning time periods, 

T on (i) is the minimum on time of generator i, 

T off (i) is the minimum down time of generator i. 

(d) Unit ramping constraints 

P(i, t) - P(i, t - 1) < UR(i) as unit i ramps up (17) 
P(i,t-l)-P(i,t)<DR(i) as unit i ramps down (18) 

where 

UR(i) is ramp up rate limit of generator i, 



DR(i) is ramp down rate limit of generator i. 

This model described above for the price-based unit 
commitment is difficult to solve. But if 'I (i, t)' is given, this 
problem became a quadratic programming problem that is easy 
to solve. Based on this observation, this paper used a greedy 
algorithm creating the units' status variables T (i, t)', and a 
quadratic program is used to maximize the profit based on 
current 'I (i, t)'. The next section will introduce the greedy 
algorithm, quadratic program, as well as how to apply them to 
solve the price-based unit commitment problem. 

IV. THE GREEDY ALGORITHM AND QUADRATIC PROGRAM 

The idea of the greedy algorithm is simple. Given a set Q l , 
the next element chosen is one that gives the greatest 
immediate increase in value, provided that such an element 
exists. Moreover, once an element is chosen, it is kept 
throughout the algorithm. Let v (Q) be a real-valued function 
defined on all subsets of N= { 1 . . . n } and consider the problem 
max {v(Q):Q c N}.The followings demonstrate the iterative 
process of this algorithm [18]. 

Initialization: Q° = § , t=l. 
Iteration t: 

Step 1: Let j t = arg max jeNVQ t-i v(Q t_1 u {j}) with ties 
broken arbitrarily. 

Step 2: If v(Q t_1 u{j t })< v(Q t_1 ) , stop. Q t_1 is a greedy 
solution. 

Step 3: If v(Q M u {j t }) > v^ 1 ), set Q 1 = Q t_1 u {j t } . 

Step 4: If Q l = N , stop. N is a greedy solution. Otherwise let 
t=t+l. 

Where v(Q) is a real-valued function defined on all subsets 
ofN={l, ,n). 

In order to specify how to apply the greedy algorithm 
cooperating with quadratic program to solve the price-based 
unit commitment problem, Table III shows an example of the 
5-unit system. Fig. 7 shows the flow diagram of the new 
algorithm for the price-based unit commitment problem The 
description of this algorithm is given below. 

TABLE III PARAMETERS FOR THE 5-UNIT SYSTEM 



Uni 


Pmax 


Pmin 


a(i) 


b(i) 


c(i) 


M 


Mi 


Ini 


Ra 


Ra 


ST 


t 












n 
on 


n 

off 


T 


mp 

Up 


mp 
Dn 


MB 
tu 


1 


12 


2.4 


24.3891 


25.5472 


0.02533 


1 


-1 


-1 


12 


12 





2 


100 


25 


217.895 
2 


18 


0.00623 


4 


-2 


-3 


50 


50 


100 


3 


260 


50 


636.46 


3.8260 


0.0089 


10 


-1 


-1 


15 


15 


500 


4 


440 


160 


669.12 


7.9215 


1.7E-06 


10 


-1 


-1 


25 


25 


603 


5 


250 


130 


590.72 


2.4435 


0.00906 


10 


-1 


-1 


12 


12 


500 



P max : Maximum Capacity (MW) 
P min : Minimum Capacity (MW) 

a(i), b(i), c(i): Fuel Consumption Coefficient (MBtu, MBtu/MW, MBtu/MW 2 ) 
Min on: Minimum on Time (Hour) 
Min off: Minimum off Time (Hour) 
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IniT: Initial Operational Time (Hour, Minus Means off) 
Ramp Up: Ramp-up Rate (MW/min) 
Ramp Dn: Ramp-down Rate (MW/min) 
St Mbtu: Start-up Mbtu 

Assign an initial state that has maximum 
on time periods for each unit 

r 

According to the initial state (I(i,t) is known), solve the 
proposed model. Assign current objective value objval 
as best objective value Bcstobj, current unit status I(i,t) 
as best unit status Bcstl(i.t), and current generation P(i,t) 
as BestP(i,t). 

4 ~ ' 

Assign contimuc=0, and assign 
Bestl(i,t) toIl(i,t) 

I 

For unit i, assign Ilto I. 




A greedy solution is found. 



Fig. 7 Flow diagram of the new algorithm for the price-based unit Commitment 

problem 

A. Initialization 

In applying the greedy algorithm to the price-based unit 
commitment, the initial set Q is assigned as all units are on 
unless it is infeasible. If it is infeasible, it must violate the 
minimum off time during some beginning time periods. The 
unit can obtain an initial maximum on time periods by turning 
off these infeasible time periods. From Table III, it can be seen 
that initial states of all units have been off for a time interval 
larger than minimum off time, therefore, all units are on for 
every time period is a feasible solution. The initial state of I (i, t) 
is shown in Table IV. 



TABLE IV THE INITIAL STATE OF I(I,T) 





t=0 




2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


21 


22 


23 


i=l 


1 
















1 




1 


1 


1 


1 






1 


1 




1 










2 


1 
















1 




1 


1 


1 


1 






1 


1 




1 










3 


1 
















1 




1 


1 


1 


1 






1 


1 




1 










4 


1 
















1 




1 


1 


1 


1 






1 


1 




1 










5 


1 
















1 




1 


1 


1 


1 






1 


1 




1 
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B. Calculate the Initial Solution 

Once an initial feasible T (i, t)' is given, this model (eq. (8) 
~ (18)) became a quadratic program. Using quadratic 
programming technique, it can determine the objective value 
'objval', and the bidding generation quantity 'P (i, t)' for each 
unit at every period. Regard the current objective value 'objval', 
current status of units T(i,t)', and generation of all units 'P(i,t)' 
as 'Bestobj', 'Bestl(i,t)' and 'BestP(i,t)' respectively. 

C. Iterative Process 

The iterative process is to let a time period of one unit off 
one time if it is on at this time period. If this change violates 
the minimum on time constraint, this is an infeasible state. But 
if it violates the minimum off time constraint, it can let it 
became a feasible state by turning off some more next time 
periods. For these given states, determine the objective value 
'objval' and 'P (i, t)'. If the current objective value 'objval' is 
greater than 'Bestobj', replace the 'Bestobj' and 'BestPri, t)' 
with current 'objval' and 'P(i,t)' respectively, and set index 
'continue' equals to 1. The index 'continue' equaling to 1 
means that the objective value is still can be improved further, 
and the iterative process should continue. Units' statuses 'I (i, 
t)' of the next iteration is based on the 'Bestlt (i, t)' of the 
current iteration. Repeating this process, until no unit turning 
off at any time period can improve the objective value. 

1 ) Quadratic program [19]: 

For a given state of units' status 'I (i, t)', the proposed 
model, becomes a quadratic program. The general quadratic 
program can be expressed as (19). 

min — x T Qx + x T c 

\ (19) 

S.T. a- x = bj i e E 
a^x < b; i e 1 

E and I are index sets for equality and inequality constraints. 
The matrix Q is symmetric and positive semi-definite (if not 
actually positive definite). 

2) Equality Constraints: 

A quadratic program is greatly simplified, and can be 
solved in closed form, if it contains equality constraints only. 
Thus we consider it in some detail the quadratic program 

min — x T Qx + x T c (20) 

2 

S.T. Ax = b 

The Lagrange necessary conditions for this problem are 

Qx + A T X + c = (21) 
Ax-b=0 

If the matrix Q is actually positive definite, then an explicit 
formula for the solution of the system can be easily derived as 
follows: From the first equation in (21) we have 

x = -Q _1 A T X-Q _1 c. 

Substitution of this into the second equation then yields 

-AQ 1 A T >.-AQ- 1 c-b = 0, 

from which we immediately obtain 

X = -(AQ" 1 A T )" 1 [AQ- 1 c + bJ (22) 

and 

x = Q- 1 A T (AQ- 1 A T r 1 [AQ-'c + bJ-Q- 1 c (23) 
= -QT 1 [i - A^AQ-'aV AQ _1 ]c + Q" 1 A T (AQ" 1 A T r 1 b. 
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3) Active Set Method 
The general quadratic program with inequality constraints 
(19) is always solved by an active set method. There is an 
especially simple version for the case where Q is positive 
definite. In this version, at iteration k a point x k is given that is 
feasible for all constraints and satisfies all the equality 
constraints of the current working set W k . The working set 
always includes the equality constraints E and possibly some of 
the inequality constraints L The quadratic program 
corresponding to the working set is then defined, by translating 
to the point x k , in the form 

1 x 

min — d k Qd k - 



TABLE VE OPTIMAL UNITS' POWER GENERATIONS (5-UNIT SYSTEM) 



Skdk 



(24) 



S.T. 



0, 



where g k =c + Qx k . 

This program has only equality constraints and can be 
solved for d k by the use of formula (23) above or by other 
numerically efficient methods. If d k =0, the current point x k is 
optimal with respect to the working set. If d k ?t 0, and x k +d k is 
feasible for all constraints, then x k +d k becomes the new x k+i . If 
x k +d k is not feasible, a search of the form x k+1 = x k +a k d k is 
made, and a k is selected as large as possible to maintain 
feasibility. The general move is therefore x k+1 = x k + a k d k , 
where 

(25) 



— minimum^ 1, 



a'd. 



V. 



TEST RESULTS 



In this paper, the 5-unit system and 36-unit system are used 
to test the proposed model. Units' parameters of the 5-unit 
system are shown in Table III. Parameters of the 36-unit 
system are available in [20]. The fuel price is 3$/MBtu. In 
order to consider uncertainties of energy prices, possible 
energy price scenarios as shown in Fig. 6 have been taken into 
account. Some test results are demonstrated below. 

In Table V, Unit 2 is shut down in the 4 th and 5 th time 
period in the first iteration for the 5-unit system. The second 
iteration determines Unit 2 should still be off at the 6 th time 
period. Executing the iteration process shown in Fig. 7, a 
converging solution will be obtained when no unit turning off 
could improve the objective function. Tables VI and VII show 
the final units' statuses and power generations respectively. Fig. 
8 shows the improving process of GENCO's profit. The final 
profit that GENCO could obtain is 503,705 dollars. 

TABLE V PARTS OF ITERATING PROCESS (5-UNIT SYSTEM) 



Iteration 
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2 


3 


4 


5 


6 


7 


8 


9 


10 


11 




Unit 




2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 




Off period 




4,5 


6 


7 


22,23 


8 


9 


20,21 


19 


10 


18 


11 





TABLE VI OPTIMAL UNITS' STATUSES (5-UNIT SYSTEM) 
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Fig. 8 the improving process of profit (5-unit system) 

Table VIII shows parts of the iterating process of the 36- 
unit system. Due to the minimum off time constraint, Unit 25 
and 26 must remain the off state till the 8 th time period initially. 
Unit 24 is shut down between the 6 th and 9 th time period in the 
first iteration. The second iteration determines Unit 23 should 
also be shut down at the same time intervals. Executing the 
iterating process shown in Fig. 7, a converging solution will be 
obtained when no unit turning off could improve the objective 
function. Table LX shows parts of the final results. Fig. 9 shows 
the improving process of the GENCO's profit. The final profit 

that GENCO could obtain is 1 .44874 x 10 6 dollars. 

TABLE VI PARTS OF ITERATING PROCESS (36-UNIT SYSTEM) 



Iteration 
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TABLE IX OPTIMAL UNITS' POWER GENERATIONS (36-UNIT SYSTEM) 
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Fig. 9 The improving process of profit (36-unit system) 

In order to verify the proposed model to obtain more profits 
than those that the crisp model can do, the data of load demand 
and energy prices between 2005/06/01 and 2005/06/07 are 
reserved. Using these real energy prices, the real profit can be 
compared between the crisp model and the proposed fuzzy 
model. Tables X and XI compare the real profits obtained by 
the regression model and the proposed fuzzy model for the 5- 
unit system and 36-unit system respectively. From these tables, 
it can be seen that the proposed model considering 
uncertainties can obtain more profits. 



TABLE X COMPARISON OF THE REAL PROFIT BETWEEN THE REGRESSION 
MODEL AND THE PROPOSED FUZZY MODEL (5-UNIT SYSTEM) 



Date 


Regression model 


The proposed model 


2005/06/01 


4.9760x10 s 


4.9805x10 s 


2005/06/02 


4.5707xl0 5 


4.5808x10 s 


2005/06/03 


4.4780x10 s 


4.4623xl0 5 



2005/06/04 


5.6284xl0 5 


5.6284xl0 5 


2005/06/05 


6.1347xl0 5 


6.1453x10 s 


2005/06/06 


8.3116x10 s 


8.3195x10 s 


2005/06/07 


9.2643x10 s 


9.3117x10 s 


TABLE XI COMPARISON OF THE REAL PROFIT BETWEEN THE REGRESSION 
MODEL AND THE PROPOSED FUZZY MODEL (36-UNIT SYSTEM) 


Date 


Regression model 


The proposed model 


2005/06/01 


1.4580x10" 


1.4729x10" 


2005/06/02 


1.2191x10" 


1.2352x10" 


2005/06/03 


1.1858x10" 


1.1878x10" 


2005/06/04 


1.6828x10" 


1.6828x10" 


2005/06/05 


1.9120x10" 


1.9408x10" 


2005/06/06 


2.9071x10" 


2.9223x10" 


2005/06/07 


3.3465x10" 


3.3733x10" 



VI. CONCLUSIONS 

In the deregulated power market, GENCOs regard 
obtaining maximum benefit as their object instead of satisfying 
the customer's demand. In this new environment, GENCOs 
should face uncertainties exiting in power market to pursue 
their profits. This paper proposes a new model considering 
uncertainties of energy prices. An efficient algorithm that 
combines the greedy algorithm and quadratic program is 
introduced to solve this model. Test results show this proposed 
model can obtain more profits than those that the crisp model 
can do. 
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